Nonclassical microwave radiation from the dynamical Casimir effect 
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We investigate quantum correlations in microwave radiation produced by the dynamical Casimir 
effect in a superconducting waveguide terminated and modulated by a superconducting quantum 
interference device. We apply nonclassicality tests and evaluate the entanglement for the predicted 
field states. For realistic circuit parameters, including thermal background noise, the results indicate 
that the produced radiation can be strictly nonclassical and can have a measurable amount of 
intermode entanglement. If measured experimentally, these nonclassicalilty indicators could give 
further evidence of the quantum nature of the dynamical Casimir radiation in these circuits. 
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Vacuum fluctuations are fundamental in quantum me- 
chanics, yet they have not so far played an active role 
in the rapidly advancing field of engineered quantum de- 
vices, e.g., for quantum information processing and com- 
munication. The main reason being that it has been no- 
tably difficult to observe dynamical consequences of the 
vacuum fluctuations [1], let alone use them for applica- 
tions. The dynamical Casimir effect (DCE) [31 [3] is a 
vacuum amplification process that can produce pairs of 
photons from quantum vacuum fiuctuations by means of 
nonadiabatic changes in the mode structure of the quan- 
tum field, e.g., by a changing boundary condition [4j [5] 
or index of refraction [6l [7] . As such it could potentially 
be applied as a source of entangled microwave photons. 

For decades the DCE eluded experimental demonstra- 
tion, largely due to the challenging prerequisite of nona- 
diabatic changes in the mode structure with respect to 
the speed of light. However, using a varying boundary 
condition in a superconducting waveguide [H [9] , the ex- 
perimental observation of the DCE was recently reported 
[lOj . This experiment also demonstrated that the dynam- 
ical Casimir radiation exhibits the expected two-mode 
squeezing [51 lllH13j . which is a consequence of a non- 
classical pairwise photon-creation process. 

The microwave radiation produced by the DCE in su- 
perconducting circuits therefore has high potential of be- 
ing distinctly nonclassical. Whether the state of a quan- 
tum field is nonclassical, or if it could be produced by a 
classical process, may be demarcated by evaluating cer- 
tain carefully-designed inequalities [2] for the field ob- 
servables (nonclassicality tests). In this paper, we apply 
such nonclassicality tests to show that the microwave ra- 
diation produced by the DCE in these superconducting 
circuits can be distinctly nonclassical, even when taking 
into account the background thermal noise [T^l [TB] and 
higher-order scattering processes. Using auxiliary quan- 
tum systems as detectors [171 IE] could be an alternative 
to directly measure the field quadratures, which could 
provide further opportunities to detect nonclassical cor- 



relations, e.g., on the single photon-pair level [l9] . 

DCE in superconducting circuits. — Superconducting 
circuits are strikingly favorable for amplifying vacuum 
fiuctuations because of their inherently low dissipation, 
which allows the vacuum state to be reached, and the 
in-situ tunability of an essential circuit element, namely 
the Josephson junction (JJ). A J J is characterized by its 
Josephson energy, and by arranging two such junctions 
in a superconducting loop - a superconducting quantum 
interference device (SQUID) - an effective tunable JJ can 
be produced. The Josephson energy of the effective junc- 
tion can be tuned by the applied magnetic flux through 
the SQUID-loop. This in-situ tunability can be used to 
produce waveguide circuits with tunable boundary con- 
ditions j20:-23 , as employed in the DCE experiment in 
Ref. [TU], and tunable index of refraction [21H2S]. Tun- 
able JJs are also essential in related DCE proposals based 
on circuit QED with tunable coupling [27] . 

The electromagnetic field confined by a superconduct- 
ing waveguide, such as a coplanar or strip-line waveg- 
uide, can be described quantum mechanically in terms 
of the flux operator ^{x,t). It is related to the voltage 
operator by ^{x,t) = J dt'V{x,t'), and to the gauge- 
invariant superconducting phase operator ip — 27r<I>/$Oj 
where $o — h/2e is the magnetic flux quantum. The flux 
field in the transmission line obeys the massless, one- 
dimensional Klein-Gordon wave equation, dxx^{x,t) — 
v~^dtt^{x, t) — 0, which has independent left- and right- 
propagating components. Using this decomposition, the 
field can be written in the form 
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a(w)e^*(-'='*'^+"*) 6(cj)e-*('=-^+"*)l , (1) 



where a(w) and &(w) are the annihilation operators for 
photons with frequency a;/27r > propagating to the 
right (incoming) and left (outgoing), respectively. Here 
we have used the notation a(— cj) = a^(a;), and k^^ = 



2 



w/u is the wavenumber, v is the speed of hght in the 
waveguide, and Zq the characteristic impedance. 

Using the previously discussed flux-tunable SQUID 
termination of the waveguide, one can produce a tun- 
able boundary condition (see also Refs. [551 122]) for the 
quantum field [Eq. ([!])], 

*(0,i) + LeffW9x$(a:,i)L=o = 0, (2) 

that can be characterized by an effective length L^s {t) — 
(<I>o/27r) /{Ej{t)LQ), where Lq is the characteristic in- 
ductance per unit length of the waveguide and Ej(t) — 
E,j[^cxt{t)] is the flux-dependent effective Josephson en- 
ergy. To arrive at this boundary condition we have 
neglected the capacitance of the SQUID and assumed 
small phase fluctuations, which is justified for a large 
SQUID plasma frequency [SI [9]. For sinusoidal modula- 
tion with frequency uid/'^TT and normalized amplitude e, 
E,j{t) — + esincj^i], we obtain an effective length 
modulation amplitude SLeS = eL^ff , where L^fj = Leff(O). 
A strong modulation (corresponding to an effective ve- 
locity Woff — SLcsUJd that is a significant fraction of the 
speed of light in the waveguide v) , results in nonadiabatic 
changes in the mode structure of the quantum field, and 
the emission of photons as described by the DCE. 

The DCE can be analyzed using scattering theory that 
describes how the time-dependent boundary condition, 
or region of the waveguide with a time-dependent index 
of refraction, mixes the otherwise independent left and 
right propagating modes [30] • The superconducting cir- 
cuits considered here were analyzed using this method 
in Refs. [HI [U], where the weak- modulation regime was 
studied analytically using perturbation theory, and the 
strong-modulation regime was studied using a higher- 
order numerical method. 

In the perturbative regime, the resulting output field is 
correlated at modes with angular frequencies co and ujd — 
Lo, i.e., symmetrically around half the driving frequency. 
This intermode symmetry is emphasized when the output 
field is written for two such correlated modes: 

.SLcff 1 , , 

b± — —a± — t ^uj^u^a', (3) 

V ^ 

where we have introduced the short-hand notation 
a± — a{u!±) and b± — fe(w±), and where u!± — 
LOd/'i ± 5lo and 5lj is the symmetric detuning. In 
this perturbation calculation, the small parameter is 
SLeffy/uj-Lu^/v ~ eLeff(0)a;d/2f . Here, even if the in- 
put field is in the vacuum state, ^aj_a±^ — 0, the out- 
put field Eq. ([3| has a nonzero, symmetric photon fiux 
bj.b±^ = {SLcs/v)^uj^iL)_, i.e., the dynamical Casimir 

radiation. Furthermore, the photons in the two modes 
have bunching-like statistics, where the probability of si- 
multaneously observing one photon in each mode is equal 
to the probability of observing a photon in one of the 

modes ib^^b+b^_b^ \ w \b±b±\, i.e., they appear in pairs. 



For finite temperatures, where thermal noise is present 
in the input field, and for not so weak modulation, when 
for example 5Leffy'w+a;_/u no longer is a small parame- 
ter, it is not obvious if or to what extent the above results 
apply. In these cases there are both classical and non- 
classical contributions to the photon flux in the output 
field, and it becomes necessary to systematically compare 
the relative importance of such contributions in order to 
tell if the resulting output field remains nonclassical or 
not. In the following, we carry out such an analysis us- 
ing nonclassicality tests and by evaluating the degree of 
entanglement in the predicted output field. 

Nonclassicality tests. — The theory of nonclassicality 
tests has been well developed in quantum optics, and 
here we briefly review the important results in the nota- 
tion introduced above for superconducting waveguides. 
We consider an operator / which is defined as a func- 
tion of the creation and annihilation operators. For the 
Hermitian operator P f it can then be shown fl4J, using 
the Glauber-Sudarshan P function formalism, that any 
classical state of the field satisfies 

(: Pf ■■) > 0, (4) 

where the condition for classicality that has been used is 
that the P function must always be non-negative. The :: 
denotes normal ordering. 

For the two-mode quadrature-squeezed states that the 
DCE is known to produce, the natural definition of / is 

fe = e'^S_ -f e~'%l + - e-'%l), (5) 

where 9 is the angle that defines the principal squeez- 
ing axis. With this definition of /g, a pure two-mode 
squeezed state is known to violate the inequality Q, 
see, e.g., Ref. [13] and references therein. This choice 
of fg is also suitable from an experimental point of view, 
since ^: flfg :^ can be evaluated from experimentally- 
accessible quadrature correlations. 

We now evaluate the quantum-classical indicator 
P f ■'^ = niin^: pfg :^ for the field state produced 

by the DCE, and discuss the conditions under which this 
nonclassicality test is violated. For weak driving, using 
output field Eq. ([3|, and a thermal input field we obtain 

{■.Pgfg-) - 2« + n">) 

- 4cos20 — V^J^(l + n*h + 71*^)^(6) 

V 

where n^} — ^aj_a±^ = {exp{hLu±/kBT) — 1)^^ is the 
thermal photon flux of the input mode with frequency 
uj±. In this case, ^: f^fg :^ is minimized by taking ^ = 0, 

and it is negative if {6Lctf/v)y/uj+uj^ ^ {n*:^ + n^^)/2, or, 
equivalently, e > 2u/(L"jfa;rf)(n*^-|-n'*)/2. This indicates 
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write the inequality /g/g :^ < as 



FIG. 1; (color online) (a) The quantum-classical indicator 
(^'- fife '-^ as a function of driving amplitude e for a range of 
6 values in the interval [0, 27r] (blue), and for ^ = (red), 
which is the optimal 6 in the perturbation regime. Due to 
the thermal input field, ^: /g/e :^ > for small e. However, 

when e is sufhciently large /^/ :) < 0, which conclusively 
rules out that the field state is of classical origin, (b) The two- 
mode squeezing (T2 as a function of the dimensionless driving 
amplitude e (red), together with the right-hand side of Eq. ^ 
(black) , which defines the boundary between the classical and 
quantum regions. Parameters: ajd/27r = 10 GHz, 5u}/uJd = 
0.15, r = 50 mK. Other parameters are from Ref. |10j . 



that the field state in the form Eq. ([s]) is distinctly non- 
classical for a vacuum input field, and potentially also for 
low-temperature thermal input fields. 

To investigate whether the nonclassical characteris- 
tics of the DCE radiation remain for realistic input 
field temperatures and when the driving amplitude is in- 
creased beyond the perturbative regime, we also evalu- 
ate ^: fife '-^ by solving the scattering problem numer- 
ically. The results of this calculation are presented in 
Fig. [l|^a) , showing that for sufficiently large driving am- 
plitude ^: f^f :^ < even at typical temperatures for su- 
perconducting circuits, and including higher-order scat- 
tering processes. We therefore conclude that nonclas- 
sical characteristics of the DCE radiation can be suffi- 
ciently robust to remain important in realistic experi- 
mental situations. Evaluating ^: f^f :) from experimen- 
tally measured field quadratures therefore appears be a 
viable method to conclusively demonstrate the quantum 
statistics of the dynamical Casimir radiation. 

The nonclassicality test in terms of ui . — To further re- 
late to the experimental demonstrations of the DCE, it is 
instructive to formulate the nonclassicality test in terms 
of the two-mode squeezing (72, which was measured in 
Ref. 10 . The two-mode squeezing is defined as 02 = 
((/_/;) -(Q_Q+))/ {{{P) + {ID + (Ql) + (g^))/2). 



where I± = (?iw±Zo/87r) 



1/2 



»05 



-"I'bl ] and 



are the volt- 



1T2 > 



[2n. 



[2n_ + 1] ' 



(7) 



age quadratures. Using this expression for 172, we can 



where n± = y^±^±J is the photon flux (thermal and 
DCE) for the output mode with frequency uj± , and where 
we have taken 9 = (f> + it /A to relate CT2 and ^: f^fg 

Equation ([t]) suggests that a non-zero two-mode 
squeezing does not necessarily imply that the field is a 
strictly nonclassical state [by the criterium of Eq. ^ 
and the current definition of the operator /]. However, 
if the magnitude of the two-mode squeezing exceeds the 
right-hand side of Eq. ([7|, the field is guaranteed to be 
distinctively nonclassical (i.e., squeezed vacuum rather 
than a squeezed thermal state). Since the expectation 
values in the right-hand side of Eq. ([7| can be measured 
experimentally, this could be a practical formulation for 
the experimental evaluation of the nonclassicality test. 

Figure shows the two-mode squeezing together 

with the boundary between the classical and quantum 
regimes, as defined by Eq. ([t]). With the parameters 
used in Fig.[T] the boundary corresponds to the squeezing 
(72 ~ 0.04. Experimental measurements |10j have demon- 
strated significantly larger squeezing for the dynamical 
Casimir radiation, but at the same time the measured 
photon flux was larger than in the current calculations 
due to the presence of low-Q resonances in the transmis- 
sion line. An increased photon flux increases the value of 
the boundary in Eq. ([7| and makes the violation of the 
inequality more demanding. However, by reducing the 
driving strength to get a lower photon flux a violation 
of the nonclassicality test Eq. ([t]) should be achievable 
with an experimental setup like the one in Ref. lOj, al- 
though increased measurement time and averaging may 
be necessary to obtain sufficient sensitivity. 

Entanglement. — The two-mode squeezing and the non- 
classicality tests discussed above demonstrate that the 
DCE radiation is nonclassical. The quantum nature of 
the radiation originates from the entanglement in individ- 
ual pairs of photons. To quantify the entanglement be- 
tween two entire modes with frequencies adding up to the 
driving frequency, we evaluate the logarithmic negativ- 
ity M j3T], which is an entanglement measure for Gaus- 
sian states that is frequently used in quantum optics, 
and recently also in microwave circuits [32) and nanome- 
chanical systems |33) . The logarithmic negativity is pos- 
itive for entangled states, and it can be calculated from 
the covariance matrix V^^ = ^ {RaRp + RfjRa), where 
R^ ~ {q-,P-, q^,p+) is a vector with the quadratures as 
elements: q± — {b± +bj.)/V2 and q± — —i{h± — ) /v^. 

The covariance matrix can be evaluated both an- 
alytically and numerically, and also constructed from 
experimental quadrature measurements. The numeri- 
cally calculated covariance matrix is shown in the in- 
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FIG. 2: (color online) The logarithmic negativity A/" as a func- 
tion of the normalized modulation amplitude e. The onset of 
nonzero A/" is eo- The parameters are the same as in Fig.[l] In- 
set: The covariance matrix for e — 0.5. The diagonal quadra- 
ture correlations correspond to the vacuum fluctuations and 
the photon flux due to the DCE and of thermal origin. The 
nonzero, off-diagonal elements correspond to the two-mode 
correlations produced by the DCE. 




FIG. 3: (color online) The region of nonclassical radiation 
(blue), visualized using — (: f^f :) (left) and the logarithmic 
negativity A/" (right), as a function of the temperature T and 
the detuning Suj, for e = 0.15 and other parameters as in 
Fig. [l] Although the nonclassical region is larger for Af than 
for (Tf^f :), Af is small in the region where (: f^f :) is non- 
negative (white), and the regions where the measures violate 
classicality with a one-cr confidence are quite similar. 



set in Fig. [2] for typical parameters. Given the covari- 
ance matrix for the two selected modes, it is straight- 
forward to evaluate the logarithmic negativity, defined 
a,s Af — max[0, — log(2j/_)], where z/_ = (a/2 — (cr^ — 
4dety)i/2/2)i/2^ and a = det ^ -f det B - 2 det C, where 
the A, B, and C are 2x2 submatrices of the covariance 
matrix V = {A,C;C'^,B). 

The logarithmic negativity for the DCE (see also 
Ref. IMl) is shown in Fig. [2] At zero temperature and 
small drive amplitudes, it is proportional to the driv- 
ing amplitude Af — eL'^gUJd/v. For finite temperatures 
and small detuning 6uj, the onset of nonzero logarithmic 
negativity is at eo ~ 2v/{L^ffLJd){n*:^n^j!^y^'^, after which 
it increases with the driving amplitude. For sufficiently 
large driving amplitude, e ^ 0.06, the quantum correla- 
tions overcome the thermal noise and the two matching 
output modes are entangled. Comparing Figs. [T] and [2] 
implies that the logarithmic negativity is a stronger in- 
dicator of the nonclassicality of the field state than the 
inequality ^ with our definition of /. This is also shown 
in Fig. [3j which visualizes the nonclassical regions as a 
function of temperature and detuning, as well as the sen- 
sitivity to uncorrelated classical quadrature noise intro- 
duced in the detector (the one-cr contour line). However, 
when taking this sensitivity into consideration, the two 
measures appear to be of similar practical usefulness. 

Conclusion. — We have theoretically investigated quan- 
tum correlations in the radiation produced by the DCE in 
a superconducting waveguide by evaluating nonclassical- 
ity tests and the logarithmic negativity. These measures 
indicates that the devices used in Ref. [10], should have 
access to regimes where the produced radiation is strictly 
nonclassical. We have formulated practical inequalities 
with experimentally obtainable observables that could be 
used to directly verify the quantum nature of the mea- 



sured radiation in future DCE experiments. We also note 
that recently two-mode squeezed states have been gener- 
ated in microwave circuits using other mechanisms, for 
example parametric amplification using the nonlinear re- 
sponse [3S] or time- varying index of refraction |2S| of 
SQUID arrays and JJs [3H[35]. The nonclassicality tests 
discussed here could also be applied to analyze the ra- 
diation produced in these experiments. We believe that 
a demonstration of a nonclassicality violation in super- 
conducting circuits, or other promising systems |37H40j . 
could pave the way to the experimental exploration of 
the continuous production of entangled microwave pho- 
tons by the DCE, and possible applications thereof in, 
for example, quantum information processing |41H43j . As 
such it could become a novel practical application of mi- 
crowave quantum vacuum fluctuations. 
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